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abstract 


research  how  fitc  vurt<k.c  grawiy  waves  travelling  over  a 

tlmd  arc  nw*bticd  by  the  presence  o t  a  submerged  pair  of  vonex  singularities, 
the  undcfvumitftg  of  vonev  flows  near  a  free  surface  has  become  a  major  goal 
in  the  hydrodynamics  vootmuruty  over  the  last  decade  due  to  us  possible 
importance  in  desiemumng  the  influence  of  vortical  ship  wai.es  on  radar  images 
of  the  ocean  surface  Although  much  research  has  concentrated  on  studying 
vortices  near  a  still  free  surface,  a  theory  for  vortical  influence  on  ambient  free 
surface  w>aves  has  not  been  addressed  The  inclusion  of  ambient  surface  waves 
is  of  fundamental  importance  in  undemanding  the  ship  wake  images  observed 
under  real  ocean  conditions  The  first  steps  in  uncovering  the  subtleties  of  the 
ambient  free  surface  wave  modification  by  submerged  vortices  are  taken  by 
concentrating  on  three  questions  First,  what  is  the  important  physical  mechanism 
of  influence  of  the  vortex  on  free  surface  ambient  waves?  Second,  in 
formulating  a  typical  vortex  model  of  the  flow  involving  a  rising  pair,  what 
parameters  are  significant  in  determining  the  form  of  the  ambient  wave 
modification0  And  finally,  how  can  this  information  be  useful  in  the  practical 
problem  of  radar  imaging  of  ship  wakes?  Careful  analytic  analysis  of  the  free 
surface  boundary  conditions  will  show  that  the  induced  surface  current  generated 
by  the  vortices  provides  the  physical  mechanism  of  ambient  wave  modification. 

A  numerical  model  is  described  which  is  utilized  to  determine  the  relative 
importance  of  parameters  derived  in  the  analytic  analysis  on  this  physical 
mechanism.  In  addition,  spectral  analysis  is  performed  to  highlight  the  use  of 
simulation  results  in  evaluation  of  conditions  applicable  to  the  problem  of  radar 
remote  sensing  of  ship  wakes. 


ADMINISTRATIVE  INFORMATION 


The  work  was  sponsored  by  the  Independent  Research  (IR)  Program  funded 
by  ONT  and  administered  by  DTRC.  The  associated  work  unit  for  this  research 
was  1543-130. 


C'HAPTtK  I  •  IV1MOIR4  IH)S 

Research  i  merest  tn  dir  ticlu^  or  » uftr  s  *o*»v ture*  t»c*i  *  ftre  vurtauc  ha> 
'ecently  been  revived  by  the  drtrxuun  of  di/i  ivg»u<i>  m  *n*gc»  ol  the  sc* 

surface  The  significance  of  these  long  rutfrow  daiV  cuewks  is  their  oouh  idcmr 
with  the  center line  '*<*'ce  region*  ui  ship*  travelling  un  the  tea  utrtiue  at  the  time 
of  radar  expos urt  foe  this  reason  they  are  gcrvrrails  referred  to  as  ia*i 
centerlines  Aft  example  erf  the  dark  centerline  irfwersTd  prsmarvly  on  jymhciK 
aperture  radar  <SAK>.  u  shown  m  fig  ta  T  e  ship  m  this  example  »s  located  at 
the  tnter section  of  the  angled  keivw  wake  lines  as  shown  diagrammat nails  in  I  ig 
lb  Other  examples  and  details  an  the  radar  operation  and  image  processing  can 
be  found  tn  t.ydcr  ct  al  *  and  hrltrer  rt  a!  ’  Although  the  cause  of  the  dark 
centerlines  is  not  well  understood,  there  is  widespread  agreement  that  the  lower 
enerf>  radar  return  in  this  region  u  a  result  of  the  climmation  of  surface  wavr* 
normally  responsible  for  construe"  jvc  interference  in  reflection  of  the  radar  signal 
Constructive  interference  occurs  when  the  component  of  the  surface  wavenumber 
vector  tn  the  look"  direction  of  the  radar  satisfies  the  following  expression 
(Wright 

k,  =  2  k  sw(8l 

where:  t»  =  2k/K.  =  surface  wavenumber 

k.  =  2n/k,  =  radar  wavenumber 
X  =  wavelength 

0  =  vertical  incidence  angle  of  the  radar  (normal  to  the  free 
surface:  0  =  0°) 
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These  waves,  termed  Bragg  Scattering  Waves,  are  typically  wind  generated,  and 
form  the  background  signal  return  level  in  SAR  images.  Dark  centerlines  in  the 
radar  images  of  moving  ships  at  sea  suggest  the  reduction  or  elimination  of  these 
Bragg  Waves  in  the  ship  wake.  This  process  has  been  suggested  by  some  to  be 
caused  by,  among  other  things,  the  turbulent  ship  wake  generated  by  the  ship. 
Measurements  by  Lindenmuth  and  Fry®  and  calculations  by  Griffin  et  al.6  of  the 
nran  and  rms  velocity  components  behind  a  1/20*  scale  twin  screw  ship  model 
indicate  that  the  mean  turbulent  energy  in  the  ship’s  wake  dies  out  over  much 
shorter  distances  than  the  obse.ved  dark  centerline.  This  has  lead  to  the  belief 
that  the  ship  generated  turbulence,  or  mean  wake  currents,  play  only  a  minor  role 
in  the  production  of  the  dark  centerline.  Observations  by  Fish  and  Blanton4 
suggest  that  the  turbulent  energy  in  the  wake  may  be  contained  in  large  scale 
vortical  structures  which  are  widely  dispersed  far  downstream  in  the  ship  wake. 
Under  these  conditions,  the  highly  intermittent  concentrations  of  vcrticity  may 
therefore  not  be  accounted  for  in  the  previous  time  averaged  measurements  and 
computations.  The  question  of  the  importance  of  vortical  wakes  on  the  dark 
surface  centerline  consequently  remains  unanswered.  Other  theories  explaining 
the  damping  of  ambient  waves  involve  the  redistribution  of  surfactants  by  mean 
vortex  induced  currents  in  the  ship  wake.  The  redistribution  process  is  therefore 
governed  by  turbulent  vortical  wake  influences  on  the  surface.  Regardless  of 
whether  this  vorticity  is  generated  by  bow  and/or  stem  wave  breaking,  propeller 
wash,  or  hull  generated  vorticity,  the  process  occurs  in  the  presence  of  ambient 
surface  waves.  These  issues  indicate  a  need  for  understanding  the  influence  of 
submerged  vortices  on  ambient  surface  waves. 


2 


In  order  to  comprehend  how  the  ship’s  wake  interacts  with  the  surface, 
researchers  have  constructed  simple  wake  models.  The  most  common  vortical 
wake  model  used  in  examining  the  interaction  with  the  free  surface  is  a  counter¬ 
rotating  vortex  pair.  This  model  is  illustrated  in  Fig  2.  The  choice  of  this 
model  is  two-fold.  First,  the  vortex  pair  may  be  used  as  a  crude  model  of  the 
wake  of  a  ship  with  either  strong  bilge  vortices  or  twin  outboard  rotating 
propellers.  Each  vortex  in  the  latter  case  represents  the  swirling  flow  generated 
by  a  propeller.  These  vortices  are  allowed  to  approach  the  surface  with  their 
mutually  induced  velocity.  The  vortex  pair  is  also  one  of  the  simplest 
configurations  for  study  and  can  be  utilized  as  a  building  block  towards 
understanding  the  much  more  complicated  turbulent  flow  interaction  wrh  the  free 
surface. 

Prior  studies  of  vortex-free  surface  interactions  focused  on  the  dynamics  of 
an  initially  flat  free  surface  under  the  influence  of  a  relatively  strong  vortex  pair. 
This  problem  is  not  directly  relevant  to  the  SAR  imaging  problem,  however,  due 
to  the  neglect  of  ambient  waves  responsible  for  the  scattering  of  radar  radiation. 
What  is  important  is  the  modification  of  the  ambient  Bragg  Scattering  Waves  by 
vortices  of  relatively  weak  strength.  The  terms  "strong"  and  "weak"  here  refer  to 
the  vortex/(free  surface)  strengths  characterized  by  a  nondimensional  parameter 
refered  to  as  the  Froude  number.  Differentiating  between  high  and  low  values  of 
the  Froude  number  determines  whether  linearized  boundary  conditions  can  be 
used  to  model  the  free  surface  evolution.  Because  weak  vortices  generate  smad 
vertical  surface  deflections,  and  no  wave  breaking  phenomena,  it  will  be  shown 
that  linearization  can  maintain  accuracy  while  providing  greater  insight  into  the 
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physics  of  the  flow.  Linearization  of  the  free  surface  also  allows  simulation  of 
simple  monochromatic  ambient  waves  without  introducing  second  order  self 
imposed  profile  modification.  The  alteration  of  these  simple  waveforms  by 
submerged  vortices  can  therefore  only  be  identified  under  linearized  conditions. 
Because  Bragg  wa  for  typical  radar  wavelengths  are  much  shorter  than  the 
depth  of  the  modelled  wake  vortices,  their  influence  on  the  vortex  paths  is 
negligible  and  will  be  ignored  in  this  study.  The  resulting  paths  of  the  vortex 
pair  will  correspond  to  the  paths  for  a  still  free  surface.  The  analytic  derivation 
of  an  alternate  form  of  linearization  under  the  above  stated  conditions  gives 
insight  into  the  dominant  physical  mechanism  responsible  for  ambient  wave 
modification  in  regions  near  submerged  free  vortices.  Free  in  this  context 
pertains  to  conditions  of  self  induced  motions  of  the  vortices.  The  motivation  for 
use  of  linearized  theory  is  based  primarily  on  the  ability  to  simulate  simple 
sinusoidal  ambient  waves  in  the  final  sections  of  the  research.  The  justification 
for  use  of  linearized  theory  will  be  given  in  chapters  2  and  3. 

Simulations  of  the  two  dimensional  vortex  pair  near  a  free  surface  with 
monochromatic  ambient  waves  are  used  here  to  determine  which  parameters  of 
the  governing  equations  are  dominant  It  will  be  shown  that  the  Froude  number 
of  the  vortices  and  their  initial  positions  play  the  most  important  roles  in 
inducing  surface  currents.  As  the  Froude  number  is  increased,  the  vortex  induced 
surface  currents  become  stronger,  resulting  in  greater  compression  of  incoming 
ambient  waves.  The  separation  distance  of  the  vortex  pair  manifests  itself  in  the 
size  of  the  disturbed  surface  region.  Other  parameters  such  as  the  ambient  wave 
steepness  and  relative  phase  resulting  from  the  analytic  study  as  independent 
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quantities,  play  secondary  roles,  as  might  be  expected  by  a  linearized 
examination. 

In  addition,  spectral  analysis  will  be  used  in  a  form  which  extracts 
characteristics  in  the  surface  similar  to  the  radar  imaging  process.  This  analysis 
will  show  that  submerged  vortices  can  give  images  with  darker  center  regions 
where  the  ambient  waves  undergo  the  modification.  This  analysis  will  be 
performed  on  a  vortex  pair  with  parameters  derived  from  typical  ship  operating 
characteristics. 
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CHAPTER  2  -  BACKGROUND  RESEARCH 

Experimental  realizations  of  the  rising  vortex  pair  model  with  a  free  surface 
were  first  carried  out  by  Sarpkaya  and  Henderson14  and  Sarpkaya15  using  a 
submerged  lifting  surface.  These  studies  showed  two  types  of  surface  signature 
of  the  trailing  vortices  which  were  termed  "scars"  and  "striations."  "Scars"  in 
this  sense  refer  to  longitudinal  depressions  in  the  free  surface  roughly  aligned  and 
positioned  just  outboard  of  each  vortex.  Scars  begin  appearing  when  the  vortices 
are  within  one  initial  separation  distance  from  the  free  surface.  "Striations"  are 
transversely  aligned  surface  deformations  which  appear  when  the  vortex  pair  are 
roughly  one  half  their  initial  separation  distance  from  the  surface.  Willmarth 
et-al.22  have  recently  carried  out  experiments  using  flaps  to  generate  two 
dimensional  (2D)  vortex  pairs  and  observe  their  impinging  on  an  initially  still 
free  surface.  Two  dimensional  computations  (also  in  Willmarth  et-al.22)  have 
been  carried  out  showing  good  agreement  with  the  experiments  until  the  point  of 
transition  of  the  vortices  to  a  turbulent  state.  Figure  3  illustrates  this  2D  vortex 
pair  configuration. 

Several  numerical  studies  of  this  2D  vortex  pair  near  a  free  surface  have 
been  carried  out  by  Marcus10,  Sarpkaya  et  al.16,  Telste20,  and  Tryggvasson19.  All 
of  these  studies  were  conducted  using  nonlinear  free  surface  boundary  conditions 
and  potential  flow.  The  free  surface  is  still  and  flat  when  the  vortices  are 
initialized.  These  studies  show  how  vortex  pairs  with  high  strength  rise  towards 
the  surface  under  self  induced  velocities  and  generate  an  upward  bulge  in  the 
free  surface  in  the  central  region  with  depressions  on  the  outer  edges.  Under 
conditions  of  very  high  vortex  strength,  they  have  shown  that  the  vortex  pair 
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may  actually  me  through  the  uuiul  surface  kxauoit  pushing  a  plume  of  fluid  up 
through  the  center  TcUar*  also  vhem-v  that  for  until  Froude  numbers  (F  = 
r7(gS*r2),  the  vortex  path  deviates  very  link  from  the  ngvd  surface  condition, 
and  that  the  surface  disturbances  are  both  small.  and  localized  to  the  region  of 
the  vortex  Trygvaxton  (1983)  ha»  aivo  simulated  the  nonlinear  free  surface 
deformation  for  a  pair  of  deformable  vortex  blobs  Blobs  are  models  of  a  vortex 
structure  composed  of  distributed  voracity  in  a  finite  region.  The  flow  is 
computed  using  potential  theory  however  and  the  vorticity  distribution  is 
discretized  into  vonex  Singularities  The  surface  deformations  at  low  Fronde 
numbers  are  insensitive  to  the  detailed  modelling  of  the  vortex  provided  the  blob 
size  is  small  compared  to  the  depth  of  the  vortex. 

For  cases  in  which  the  vortex  disturbance  to  the  surface  is  weak  in  some 
sense,  Wehausen  and  Lju tone21  give  an  exact  solution  for  a  linearized  free  surface 
following  the  classical  assumptions  used  for  surface  wave  theory.  Both 
Marcus10  and  the  author  have  exercised  numerical  evaluation  of  this  solution  (a 
double  integral  representing  the  Green’s  Function  over  the  surface  in  time)  over 
short  time  periods.  The  author  has  derived  an  alternate  form  of  the  classical 
solution  (Wehausen  and  Laitone21),  given  in  Appendix  A  of  Fish1,  particularly 
suited  for  conditions  of  shallow  vortices.  However,  both  classical  formulations, 
have  been  found  by  the  author  to  be  inadequate  under  conditions  of  self-induced 
motion  of  the  vortices.  This  deficiency  lies  in  an  improper  zero"'  order  condition 
around  which  the  linearization  is  performed,  and  is  described  in  further  detail  in 
the  next  section. 
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Recent  viscous  flow  simulations  of  the  vortex  pair  approaching  an  initially 
still  free  surface  have  been  carried  out  for  laminar  conditions  including  nonlinear 
free  surface  boundary  conditions  by  Ohring  and  Lugt12.  The  cunent  limitation  of 
low  Reynolds  number  of  these  calculations  prevent  them  however  from 
simulating  flows  of  the  typical  ship  scale.  In  addition,  inclusion  of  surface 
waves,  with  their  intrinsically  large  Reynolds  numbers,  still  appears  to  be  only  a 
long  range  goal  in  their  work. 

Studies  of  wave  modification  by  steady  currents  have  also  been  pursued  in 
recent  years.  Important  differences  exist  however  between  these  works  and  the 
present  study.  Cooper2  and  Griffin  et  al.6  used  conservation  of  wave  action  to 
examine  the  modification  of  short  waves  by  surface  current  patterns.  The  surface 
current  model  used  consisted  of  an  axial  velocity  defect  simulating  the  wake  of  a 
ship.  Transverse  currents  were  assumed  negligible  as  were  vertical  current 
gradients.  The  orientation  of  the  vortices  as  shown  in  Fig  2  generate  transverse 
currents,  and  allow  vertical  current  gradients  as  defined  by  the  flowfield  of  the 
vortex  pair.  The  most  important  difference,  however,  is  the  allowance  of 
significant  transverse  current  gradients  in  the  direction  of  wave  propagation.  As 
will  be  found  later,  this  condition,  in  conjunction  with  the  short  time  scales  of 
flow  evolution  for  moderate  vortex  strengths,  prevents  the  use  of  wave  action 
conservation  theory.  Surface  wave  modification  by  vertical  current  gradients  have 
been  studied  by  Simmen  and  Saffman17  and  by  Teles  Da  Silva  and  PeregrineI, 
(1988).  These  studies  concentrated  on  steady  state  alterations  of  deep  and 
shallow  water  waves,  respectively,  by  a  linear  vertical  current  gradient.  A 
particularly  important  result  of  these  works  is  the  finding  that  the  kinetic  and 
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potential  energies  of  the  modified  wave  are  unequal,  suggesting  that  wave  energy 
spectra  calculated  from  surface  height  measurements  may  be  inaccurate  when 
current  gradients  are  present.  Recent  measurements  by  Fry  and  Huang5  of  the 
modification  of  surface  waves  by  submerged  turbulence  and  a  rising  (somewhat 
two-dimensional)  jet  showed  that  modification  was  present,  but  give  no  details. 

In  the  present  work,  the  transient  response  of  the  ambient  waves  will  be 
examined  rather  than  the  steady  state  condition. 

In  order  to  address  the  specific  flow  conditions  observed  behind  real  ships, 
this  research  takes  the  first  step  in  understanding  the  flow  interaction  between 
submerged  voracity  and  an  ambient  surface  wave  field.  The  influence  of 
submerged  vortex  singularities  on  ambient  surface  waves  will  be  characterized  by 
numerically  simulating  the  transient  flow  of  an  impulsively  started  vortex  pair 
below  an  ambient  wave  field  Figure  4  shows  a  diagram  of  the  flow  being 
examined.  Both  experimental  observations  and  numerical  calculations  indicate 
that  vortex  pairs  of  practical  interest  in  the  influence  of  the  dark  radar  return 
regions  of  the  wake  are  of  relatively  low  Froude  number.  One  expects  therefore 
that  a  linearized  free  surface  boundary  condition  can  be  used  to  sustain  a  simple 
monochromatic  ambient  wave  for  modification.  The  classical  linearizations  used 
by  previous  researchers  however  contain  no  analytical  means  of  interaction 
between  the  vortices  and  the  ambient  waves.  The  solutions  are  simple 
superpositions  of  the  ambient  wave  profile  with  the  surface  solutions 
corresponding  to  the  initially  still  free  surface.  A  new  method  of  linearization 
must  be  developed  to  include  interaction  terms  in  the  boundary  conditions 
between  the  vortices  and  the  ambient  waves.  This  method  is  derived  by 
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examining  the  physical  mechanism  of  the  interaction,  and  developing  the  propei 
choice  of  zero*  order  solution  around  which  to  linearize.  The  resulting 
interaction  theory  is  described  in  the  next  section. 
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Figure  3  -  Typical  Vortex/Free  Surface  Problem 


Figure  4  -  Ambient  Wave  Modification  by  Submerged  Vortices 


CHAPTER  3  -  PHYSICAL  MECHANISM  FROM  NEW  LINEARIZATION 

The  development  of  a  theory  illuminating  the  physical  mechanism 
responsible  for  ambient  surface  wave  modification  starts  with  examination  of  the 
flow  equations.  Free  surface  flows  with  waves  arc  typically  of  such  high 
Reynolds  number  that  the  flow  can  be  considered  mviscid  (Newman51)  The 
induced  velocities  in  the  vortex/free  surface  flows  of  concern  here  arc  sufficiently 
small  that  compressibility  effects  can  also  be  ignored  (Newman11).  A  two 
dimensional  (2D)  domain  will  be  used  here  primarily  for  its  simplicity  in 
handling  the  quasi-2D  vortex  pair  ship  wake  model  described  in  the  introduction 
and  illustrated  in  Fig  2.  The  quasi- 2D  modelling  assumption  is  that  the  velocity 
gradients  along  the  axis  of  the  ship  are  small  and  that  a  2D  flow  realization  in  a 
plane  normal  to  the  ship  axis  may  be  integrated  in  time  and  stacked  in  space 
using  the  constant  ship  speed  to  represent  the  three  dimensional  problem. 

The  ambient  waves  arc  propagated  from  the  side  of  the  2D  domain  as 
shown  in  Fig  4.  These  waves  represent  ambient  waves  moving  transverse  to  the 
ship  track.  Because  of  the  quasi-2D  nature  of  the  ship  wake  model,  modelling 
waves  propagating  in  directions  much  different  from  this  lateral  condition  should 
not  be  assumed  practical. 

The  free  surface  boundary  conditions  are  being  linearized  for  two  reasons. 
The  first  reason  for  linearizing  the  free  surface  is  to  allow  stable  solutions  for  the 
propagation  of  sinusoidal  ambient  waves  in  the  absence  of  the  vortices.  Vortex 
induced  modifications  to  this  simple  ambient  wave  form  may  then  be  easily 
detected  and  quantified  using  spectral  analysis.  The  second  reason  for  linearizing 
the  free  surface  boundary  conditions  is  to  obtain  simpler  equations  to  solve.  If 
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the  linearization  is  valid,  the  equations  will  reflect  the  correct  physics  for  the 
problem  and  will  be  easier  to  interpret. 

This  section  will  develop  a  new  form  of  linearized  boundary  conditions 
using  a  more  appropriate  form  of  zero'1'  order  solution  than  that  used  in  classical 
free  surface  flow  linearizations.  This  new  linearization  method,  which  will  be 
refered  to  as  the  alternate  linearized  method  (ALM),  properly  accounts  for  the 
interaction  of  submerged  vortices  moving  under  self  induction  with  the  free 
surface.  Once  the  ALM  is  developed,  comparisons  are  made  with  both  the 
classical  linear  and  the  fully  nonlinear  theories  of  vortex  interaction  with  a  free 
surface.  Additional  insight  and  justification  for  the  ALM  will  be  given  in  this 
comparison.  The  physical  mechanism  responsible  for  the  modification  of  ambient 
surface  waves  also  becomes  apparent  in  this  comparative  analysis. 

All  of  the  equations  used  in  this  study  are  nondimensionalized.  Subscripted 
notation  will  also  be  used  to  signify  the  partial  differential  operator  unless 
otherwise  noted. 

For  example:  =  d§/dy 

3.1  BOUNDARY  CONDITION  DEVELOPMENT 

Under  the  above  described  conditions,  the  flow  can  be  represented  using  a 
velocity  potential,  the  gradient  of  which  gives  the  fluid  velocity  field.  This 
velocity  potential  must  be  analytic  in  the  fluid  domain,  and  therefore  satisfy 
Laplace’s  equation  (Newman11): 

4>«  +  4>yy  =  o  (for  two  dimensional  flow) 
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The  singular  representation  of  a  2D  point  vortex  must  therefore  be  removed 
from  the  fluid  domain  by  properly  choosing  the  boundary  of  the  flow  to  include 
a  vanishingly  small  enclosure  of  this  singularity. 

Additional  equations  describing  the  flow  are  obtained  by  enforcing  boundary 
conditions  at  the  edges  of  the  simulation  domain.  The  fully  nonlinear  boundary 
conditions  at  the  free  surface  are  given  by  (Newman11): 

Dynamic  Boundary  Condition:  (on  y  =  rj) 


<j>,2(x,T),t)  +  4>/(xJl,t)  q(x,t) 

♦.OW)  + - - - +  _  =  0 

2  F 

Kinematic  Boundary  Condition:  (on  y  =  q) 

D 

—  (y-Tj)  =  ^(x.q.t)  -  T)t(x,t)  -  0,(X,T|,t)  T),(x,t)  =  0 


(la) 


(lb) 


where:  T 

F  = - 

(g  D*r 


tT 


D  =  characteristic  length  scale  (usually  vortex  depth) 
r  =  vortex  circulation 
g  =  gravitational  acceleration 
0  =  velocity  potential 

rj  =  free  surface  elevation  above  still  surface  level 
(note:  ’  represents  dimensional  quantity) 


Vortex  Velocity  Condition: 


V  =  induced  velocity  from  other  vortex  and  images  in  free 
su’face 


These  equations,  in  addition  to  Laplace’s  equation  in  the  fluid  domain  must  be 
satisfied  for  the  flow  simulation.  The  linearization  process  is  applied  to  the  free 
surface  boundary  conditions  for  the  purpose  of  simplifying  the  mathematics  and 
numerical  flow  simulation  model  while  retaining  the  dominant  features  of  the 
flow. 
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The  first  step  in  a  linearization  procedure  is  the  choosing  of  both  the  state 
or  condition  around  which  one  wants  to  hncanzc,  and  the  small  perturbation 
parameter.  The  linearized  solution  will  therefore  be  composed  of  this  baseline 
solution  plus  a  power  series  of  solutions  m  the  small  perturbation  parameter 


4>  =  <>0  +  e4>, 

+  e^j  +  ••• 

(2ai 

ii 

+ 

3 

+  e^rij  + 

(2b) 

Each  <{>,  must  satisfy  Laplace's  equation  V;o.  =  0  in  the  fluid  domain  The 
choice  cf  baseline  solution  around  which  to  linearize  is  governed  b>  the  dewrr  to 
minimize  the  value  of  the  perturbation  senes,  hence  the  desire  to  use  the  closest 
known  approximate  solution  to  the  problem  being  addressed 

In  the  classical  linearization,  one  assumes  that  I"  is  the  small  perturbation 
parameter.  This  implies  that  the  disturbances  in  0  and  arc  small  The  bav; 
state,  or  zero*  order  solution,  therefore  is  *  tv  *  0  everywhere  The 
derivation  of  the  resulting  classical  linearized  boundary  conditions  can  be  found 
in  Appendix  B  of  Fish’  The  resulting  classical  linearized  method  (CLMj  results 
in  linearized  free  surface  boundary  condition:  given  by  (note  6  *  * 

fio+eir,) 


Dynamic: 

© 

ii 

§ 

♦ 

♦ 

(on  y  =  0) 

(3a) 

Kinematic: 

♦  n,  *  o 

(on  y  =  O' 

(3b) 

In  the  CLM  ther  are  two  distinct  governing  flow  parameters.  One  of  these,  the 
speed  of  the  vortex,  is  assumed  to  be  of  order  one;  and  the  other,  the  Froude 
number  based  on  T  is  be  assumed  to  be  of  order  e.  Further  discussion  of  these 
parameters  and  their  implications  will  be  made  in  the  section  "Comparison  of  the 
Methods." 
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This  chapter  shows  that  a  new,  more  representative  soluoon  may  be  used  as 
the  baseline,  or  0“  order  solution  to  provide  greater  insight  into  the  interaction  of 
vortices  with  a  free  surface 


3.1.1  Alternate  Linearization  .Method  (ALM) 

The  ALM  is  developed  around  a  more  meaningful  O'*  order  solution  for 
flows  with  submerged  vortex  singularities  of  low  Froodc  number  (F)  In  this 
ca*c  the  baseline  solution  is  the  resulting  flow  in  the  limit  as  F.  the  small 
perturbation  parameter  goes  to  zero  This  condition  corresponds  to  having  the 
gravitational  acceleration  so  much  greater  than  the  vortex  strength  divided  by  the 
vortex  depth,  that  the  surface  boundary  conditions  approach  those  of  a  ngid  wall 
The  zero*  order  solution  therefore  corresponds  the  ngid  wall  condition,  and  is 
formulated  as  the  sum  of  the  potential  of  the  vortex  and  its  image  above  the 
surface:  (in  dimensional  form) 

0o  =  Olid  [ln(z-c)  -  ln(z-c’)] 

Here,  z  =  x+iy  and  c  is  the  complex  location  of  the  vortex  below  the  ngid 
surface.  The  c*  represents  the  complex  conjugate  or  image  of  c  when  the  rigid 
boundary  is  located  at  y  =  0.  The  corresponding  surface  elevation  for  the  ngid 
wall  condition  is,  by  definition,  =  0.  This  can  be  confirmed  by  substitution  of 
F  =  0  into  the  Dynamic  Boundary  condition.  The  perturbation  expansion 
parameter  (e)  will  represent  the  value  of  the  Froude  number  squared  for 
simplicity,  so  this  method  can  be  assumed  valid  only  for  small  values  of  F. 
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Because  the  nonlinear  boundary  conditions  are  satisfied  on  the  deflected 
surface  elevation  position,  and  in  linearizing  one  would  like  to  satisfy  the 
boundary  condi  dons  on  y  =  0  (Tj0  =  0),  one  must  expand  the  boundary  conditions 
(1)  in  a  power  series  in  (an  order  e  quantity)  for  terms  evaluated  at  y  =  0. 

For  example,  take  a  single  term  from  the  dynamic  boundary  condition: 

4>,(x,Ti,t)  =  4>,(x,0,t)  +  T)  4>„(x,0,t)  +  Tf  $w(x,0,t)/2  +  ••• 

The  resulting  expanded  boundary  condition  equations  are  shown  below:  (keeping 
only  two  terms  in  each  Taylor  scries) 

Dynamic:  <j>.  +  fl  +  1/2  ((<>,  +  (j^)2  +  (<J>y  +  <J>„,)2)  +  fl/e  =  0  (4a) 

Kinematic:  >y  -  +  ($,  +  $„)(n,  +  A  fl„)  =  0  (4b) 
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The  <j)  and  T)  expansions  (2)  are  next  substituted  into  the  above  equations  to 
obtain  the  following  expressions: 

Dynamic: 

(<j>ot+  £01(  +  e^z)  +  (erii  +  e2^)  (0o  +  £01/  ^0^  + 

0.5  [(0o  +  £0i  +  £"02)  +  (£T1,  +  £^2)  (0o  +  £0i  +  £202)]2  + 

0.5  [(4>0  +  £<>,+  £*02)  +  (£Tli  +  £^2)  (00  +  £01  +  £202>]2  + 

y  y  y  yy  yy  n 

(EH,  +  &\2)/e  =  0 
Kinematic: 

-(0o  +  £0i  +  £"02)  -  (£Tli  +  £^2)  (0Oyy+  £01,+  e*04+  (en'«+  ^  + 

(en,  +  e2^)  (£tji^+  +  (4>0  +  £0,/  e^z)  (eni,+  £^2)  + 

(£11,  +  e2^)  (en,x+  e27]2)x((p0  +  £0lx+  £202)y+ 

(en,  +  E^z)  (<pox+  £0,x+  e^z)  (£T1,  +  e'Tlz)  + 

(en,  +  e2^)  (<j>0]  +  £0,^+  ^zHeni  +  £^4  =  0 
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These  equations  can  then  be  broken  down  into  sub- equations  in  each  power 
of  £.  The  Dynamic  Boundary  condition  is  multiplied  through  by  £  to  remove 
the  e  from  the  denominator  of  the  last  terms  on  the  left  hand  side.  The  resulting 
equations  may  then  be  compared  at  equivalent  flowers  of  £  as  follows: 


Dvnamic: 

£°: 

si 

o 

II 

O 

(5a) 

El: 

0o  +  0.5  0O2  +  T),  =  0 

t  % 

(5b) 

E2: 

01  +  4>0  01  +  ^2  =  0 

1  X  X 

(Sc) 

£3: 

02  +  01  +  0o  02  +  0.5  ({I,2  +  T|,  4>0  0,  +  0,2  +  fiy  =  0 

*  ty  *  *  l  »  ry  y 

(5d) 

Kinematic: 

ef: 

-0o  =  O 
y 

(6a) 

e1: 

-01  +  Tl,  +  00  Til  -  Tl,  00  =  o 

(6b) 

y  «  *  *  yy 


-<t>2  -  ti,  <}>,  +  t]2  +  if,  T|j  +  <j>0  r\2  +  0,  tj,  +  q,  0O  il,  -  r\2  0O  =  0  (6c) 

y  yyt  iy**»*  I  *y  yy 

The  first  order  problem  (Eq.  Sb)  shows  that  the  1“  order  free  surface 
elevation  can  be  computed  as  a  by  product  of  the  known  zerolh  order  potential. 
The  known  t),  and  0O  can  then  be  substituted  into  the  Laplace’s  equation  and 
boundary  condition  6b  to  determine  0,.  Known  0o,  0„  and  T),  may  then  be  used 
in  boundary  condition  Sc  to  determine  Tfe.  This  method  of  cycling  through  the 
boundary  condtion  equations  of  increasing  order  for  more  accurate  solutions  is  a 
characteristic  of  quasi-steady  problems.  Solving  the  unsteaddy  free  surface 
problem  by  this  method  is  extremely  combersome,  however.  The  prefered 
solution  method  is  to  solve  the  Laplace  equation  at  an  instant  in  time  and 
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integrate  this  solution  forward  in  time  using  the  boundary  conditions.  This  time 
integration  method  of  solution  is,  in  fact,  required  to  properly  incorporate 
interaction  terms  between  the  vortex  induced  velocities  and  the  ambient  waves. 

By  eliminating  all  terms  of  order  e2  or  greater,  the  following  expressions 
are  obtained: 

Dynamic: 

0o  +  e  +  $oV2  +  r|»  =  o  (7a*) 

l  t  X 

Kinematic: 

-<4>o  +  e  0,)  +  e  q,  +  e  0O  q,  -  e  q,  0O  =  0  (7b) 

y  7  t  X  X  yy 

One  can  now  obtain  a  solution  of  tj  to  second  order  by  adding  terms  from  the 
left  side  of  Eq.  5c,  who’s  sum  is  zero,  to  the  dynamic  boundary  condition  7a’. 
The  resulting  Dynamic  boundary  condition  is  given  by: 

Dynamic: 

00  +  £  4>,  +  4>oV2  +  E  0O  +  Ti,  +  e  t)2  *  0  (7a) 

I  t  X  XX 

In  order  to  compare  these  equations  to  Eq.  3  from  the  Classical  Linearization 
Method  (CLM),  they  should  be  rewritten  in  terms  of  the  variables  0  and  q 
(where:  0  =  0o+E0i  and  q  =  eq^e^,  note:  0O  is  different  in  CLM  from  ALM). 
Several  order  e2  terms  must  be  added  to  Eq.  7  for  this  rewriting.  For  example, 
0.5  0U2  is  added  to  the  dynamic  boundary  condition  7a,  and  (0Qlq2+q20oyy)  is 
added  to  the  kinematic  boundary  condition  7b.  This  adding  of  higher  order 
terms  is  legitimate  as  long  as  the  final  result  is  only  considered  accurate  to  the 
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original  order:  0  =  order(e),  *1  =  order(e2).  The  following  rewritten  equations 
result: 

Dynamic: 

0,  +  0.5  0„2  +  ltfe  =  0  (8a) 

Kinematic: 

-4>,  +  +  <t>o  T\  -  n  4>o  =  0  (8b) 

*  *  rr 

3.1.2  Comparison  of  the  Methods 

The  alternate  linearized  method  (ALM)  should  now  be  compared  with  both 
the  classical  linearization  method  (CLM)  and  the  nonlinear  method  to  determine 
the  flow  conditions  under  which  each  of  these  methods  produces  the  most 
accurate  modelling.  The  examination  starts  with  evaluation  of  the  characteristic 
parameters  of  the  flow  being  simulated.  For  a  single  vortex,  a  dimensional 
analysis  of  the  flow  indicates  that  two  independent  parameters  can  be  defined  as 
shown  below: 

T  V 

(g  D1)10  =  (g  D)m 

F  represents  the  traditional  Froude  number  based  on  the  strength  of  both 
the  vortex  and  the  gravitational  field,  and  a  characteristic  length  scale  such  as  the 
depth  of  submergence.  This  is  the  small  parameter  used  in  the  ALM.  As  F 
increases  the  vortex  exerts  a  greater  influence  on  the  free  surface  by  a 
combination  of: 

(a)  strong  vortex  circulation  strength 

(b)  shallow  submergence  depth 

(c)  weak  gravitational  acceleration  strength 
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This  parameter  must  be  small  for  validation  of  either  method,  due  to  their 
assumptions  that  the  surface  disturbance  is  a  small  quantity. 

The  parameter  G  is  an  independent  parameter  relating  the  translational 
speed  of  the  vortex,  the  gravitational  strength,  and  a  characteristic  length  scale. 
Note  that  nothing  mathematically  has  yet  been  specified  concerning  the  speed  of 
the  submerged  vortex.  The  extra  terms  of  the  ALM  boundary  conditions  should 
therefore  be  examined  for  varying  G  while  keeping  F  small. 

For  large  G,  the  physical  problem  involves  a  vortex  which  is  moving  at 
high  speed/depth  ratio.  This  is  typical  of  moving  lifting  surfaces  which  are 
forced  at  a  much  higher  velocity  than  that  induced  on  the  foil  by  the  surrounding 
fluid.  Under  these  conditions,  the  time  derivative  of  the  potential  on  the  free 
surface  may  be  much  larger  than  the  spatial  derivative.  Looking  at  Eq.  8,  the 
additional  terms  have  spatial  derivatives  of  second  order  and  can  be  assumed  to 
be  small  compared  with  the  time  derivative  and  first  order  spatial  derivative 
terms.  This  indicates  that  the  CLM  boundary  conditions  given  by  3  are  good 
approximations  for  high  values  of  G. 

If  G  is  small,  the  spatial  derivatives  may  be  of  the  same  order  as  the  time 
derivatives  and  the  additional  terms  given  by  the  ALM  boundary  conditions  may 
be  significant.  To  estimate  the  significance  of  these  extra  terms  for  a  "free" 
vortex,  the  following  analysis  is  given.  For  small  F,  the  velocity  of  the  vortex 
can  be  approximated  by  the  velocity  induced  on  its  position  by  its  image  vortex 
above  the  surface.  This  velocity  can  be  expressed  nondimensionally  by  V  = 
l/(4rcD).  In  order  to  have  the  additional  term  in  the  Dynamic  boundary 
condition  be  of  same  order  as  the  time  derivative  term  we  can  change  to  a 
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constant  velocity  coordinate  system  moving  with  the  vortex  to  express  the  time 
derivative  term  as:  4>,  =  V  <J>M.  Note  that  for  small  F,  4>,  can  be  approximated  by 
l/(2itD)  and  V  can  be  approximated  by  l/(4rtD).  The  resulting  expressions  for 
and  <b*2/2  are  identical: 

♦t  *  V  «  l/tfJtD2)  $,72  =  lA8it*D») 

The  effect  of  neglecting  the  extra  ALM  terms  in  computation  of  low  Froude 
number  surface  deflections  with  "free"  vortices  will  be  shown  in  the  next  section 
describing  the  interaction  mechanism.  Figure  9  will  show  the  disagreement 
between  the  two  methods  of  linearization  and  the  erroneous  prediction  of  an 
upward  hump  rather  than  a  downward  dip  in  the  surface  over  the  vortices. 

The  above  discussion  shows  that  under  "free"  vortex  conditions  and  low 
values  of  P/gD3,  the  extra  terms  given  by  the  ALM  contribute  significantly  to 
the  calculation  of  the  surface  deformation.  If  one  places  an  additional  constraint 
on  the  flow  by  prescribing  the  way  the  vortex  will  move,  then  the  originally 
independent  G  parameter  becomes  dependent  on  F  through  this  constraint  This 
comparison  shows  that  if  the  vortex  moves  as  a  "free"  singularity,  i.e.,  it  moves 
with  velocity  prescribed  by  the  gradient  of  the  potential  (excluding  the  local 
singularity)  evaluated  at  its  position,  then  G  is  small  when  F  is  small. 

Linearized  modelling  of  free  surface  flows  with  submerged  "free”  vorticity  should 
therefore  apply  the  ALM  instead  of  the  CLM  for  accurate  simulation  of  the  flow. 
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3.2  INTERACTION  MECHANISM  -  INDUCED  SURFACE  CURRENT 
In  order  to  evaluate  the  importance  of  the  additional  terms  present  in  the 
ALM,  a  numerical  simulation  code  was  performed  with  and  without  these  terms. 
The  simulations  involve  both  still  and  wavy  free  surfaces.  The  still  surface  cases 
are  performed  for  comparison  with  nonlinear  results  obtained  by  Telste30.  The 
wavy  surface  case  is  simply  a  sample  of  the  many  cases  explored  in  the 
parametric  study  of  the  next  section,  and  is  used  to  highlight  the  influence  of  the 
terms  representing  the  vortex  induced  surface  current  Before  describing  these 
results  however,  a  brief  description  of  the  numerical  technique,  known  as  the 
Boundary  Integral  Technique  (BIT),  is  given.  More  detailed  discussion  of  the 
BIT  is  given  in  Appendix  C  of  Fish5. 

3.2.1  Numerical  Simulation  Method:  (BIT) 

The  Boundary  Integral  Technique,  or  BIT,  is  a  commonly  used  numerical 
method  for  potential  flow  study.  The  BIT  conceptually  splits  the  solution  method 
into  two  distinct  parts.  The  first  part  is  the  solution  of  Laplace’s  equation  in  the 
fluid  domain  given  the  distribution  of  either  <J>  or  y  on  the  boundary.  The 
second  part  of  the  solution  is  the  extrapolation  of  new  values  of  <j)  or  \jr  on  the 
boundaries  using  the  unsteady  boundary  conditions. 

In  addressing  the  first  part  of  the  BIT  solution  Cauchy’s  integral  formula  is 
applied  to  the  fluid  domain.  Cauchy’s  Theorem  states  that  a  nonsingular  function 
in  a  specified  domain  is  completely  defined  by  its  values  on  the  bounding  surface 
through  Cauchy’s  integral  formula.  The  BIT  is  utilized  in  this  study  in  the 
following  manner.  First,  one  constructs  the  complex  potential  (3(z,t)  = 
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$(z,t)+i\|/(z,t)  where  z  =  x+iy.  If  P  is  analytic  then  Cauchy’s  integral  formula 
states  that 


( 

J  Q 


P 


dz  =  0 


(7) 


'o  z  -  Zo 

where:  is  the  closed  contour  integral  along  the  boundary  Q  and:  z„  is  outside 

the  boundary.  Figure  5  shows  that  the  boundary  Q  is  composed  of  both  the  free 
surface  and  all  other  boundaries.  Since  vortex  singularities  exist  inside  our 
boundary,  infinitesimally  small  circular  boundaries  must  be  added  around  them. 
The  practical  implication  of  this  is  that  the  known  complex  potential  due  to  each 
singularity  must  be  subtracted  from  the  total  complex  potential  to  obtain  an 
analytic  function.  The  resulting  perturbation  potential  is  solved  using  the  BIT 
and  added  back  to  the  singularity  potentials  to  obtain  new  total  potential  values. 

P  here  will  therefore  refer  to  the  analytic  complex  perturbation  potential. 

In  the  limit  as  Zq  approaches  Q  from  the  outside,  a  small  semi-circle  may 
be  formed  in  O  at  Zo  whose  contribution  to  7  can  be  extracted  forming  the 
following  equation: 

P 

_ dz  +  iOop(Zo,t)  =  0 

Q  Z'Zo 

where:  Ja  is  a  Cauchy  principal  value  integral  around  Q,  and  Oo  is  the  included 
angle  between  left  and  right  £2  tangents  at  z„. 

One  may  discretise  the  contour  into  N  panels  for  numerical  integration  of 
this  equation  evaluated  at  N  Zo  positions  on  the  boundary.  The  result  is  N 
equations  and  2N  unknowns  corresponding  to  the  0  and  \|/  values  at  the  discrete 
boundary  points.  If  one  can  prescribe  values  for  either  $  or  y  at  each  boundary 
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point,  the  missing  component  may  be  solved  for.  On  the  free  surface,  one  may 
prescribe  a  known  initial  distribution  of  <j>.  On  all  other  boundaries  y  will  be 
prescribed.  The  common  points  between  free  surface  and  other  boundaries  will 
therefore  have  both  0  and  V/  known.  This  results  in  N  ~  2  unknowns  and  one 
may  choose  to  solve  either  the  real  or  imaginary  part  of  Eq.  7.  Fewer  unknown 
quantities  are  necessary  if  the  real  part  is  used  on  the  Zq  points  corresponding  to 
the  free  surface  and  the  imaginary  part  is  used  on  the  points  associated  with 
the  other  boundaries.  The  forms  of  these  equations  are  given  below.  The 
resulting  N  -  2  integral  equations  are  approximated  using  trapezoidal  spatial 
summation  and  may  be  solved  for  the  N  -  2  unknowns  by  matrix  inversion,  the 
details  of  which  are  given  in  Appendix  C  of  Fish*. 

The  resulting  $  and  y  values  represent  equilibrium  solutions  to  Laplace’s 
equation  for  a  single  time  step.  These  values  are  numerically  differentiated  as 
necessary  using  three  point  central  differencing  for  substitution  into  the  free 
surface  boundary  conditions.  The  boundary  conditions  are  used  in  a  trapezoidal 
time  integration  form  as  follows: 

Kinematic:  T](t+At)  =  *n(t)  +  At  (0y  -  $0  "H,  +  "H  <J>0  ) 

*  yy 

Dynamic:  (j)(t+At)  =  4>(t)  +  At  (.5<J>„2  +  <J>0  4>i  +  T\/F) 

X  XX 

to  advance  the  solution  to  the  next  time  step.  This  time  integration  is  actually 
carried  out  in  a  predictor/corrector  scheme  for  greater  stability.  The  details  of 
this  operation  are  also  given  in  Appendix  C  of  Fish3. 

Several  special  considerations  in  simulating  free  surface  flows  deserve 
mentioning.  The  first  is  the  handling  of  radiation  conditions.  Radiation 
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conditions  typically  prevent  solutions  which  involve  wave  propagation  into  the 
flow  domain  from  outside.  Wave  damping  regions  near  the  edges  of  the  free 
surface  have  been  used  by  some  to  prevent  "reflections"  of  outgoing  waves. 
Telste20  gives  an  excellent  example  of  this.  Damping  regions  are  not  used  in  this 
research  because  of  their  detrimental  effect  on  the  ambient  waves  being 
simulated.  Ambient  wave  reflections  are  minimized  by  prescribing  boundary 
conditions  at  the  edges  of  the  free  surface  corresponding  to  the  undisturbed 
ambient  wave  potential  and  stream  function.  The  combination  of  this  prescription 
of  wave  potential  on  the  edges  with  an  initial  distribution  of  both  wave  potential 
and  elevation  at  t  =  0  provide  the  ambient  wave  energy  needed  in  the  simulation. 
The  ambient  wave  model  used  here  is  a  monochromatic  sinusoidal  deep  water 
wave.  Over  the  time  periods  simulated  here,  this  technique  proved  adequate. 

The  small  amount  of  wave  reflection  caused  by  coupling  of  vortex  potential  with 
ambient  wave  potential  on  the  edge  of  the  free  surface  was  small  and  appeared 
to  have  little  influence  on  the  main  interaction  region  over  the  time  scales 
simulated. 

The  computer  code  used  in  this  research  was  essentially  written  from 
scratch  in  FORTRAN.  The  author  used  "canned"  subroutines  for  performing  the 
LU  Decomposition  of  the  geometric  influence  coefficients  of  the  boundary  in  the 
Cauchy  integral,  and  the  numerical  differentiation  of  variables  for  the  free  surface 
boundary  condition  solution. 
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3.2.2  Numerical  Accuracy 

The  simulation  program  is  evaluated  for  numerical  error  by  simulating  the 
known  solution  for  linearized  deep  water  waves  with  extremely  small  vortex 
strength.  Zero  vortex  strength  could  not  be  simulated  by  virtue  of  the 
nondimensionalization  scheme  chosen  for  the  boundary  conditions.  The  reason 
for  this  is  that  coefficients  in  both  the  wavemaker  boundary  conditions  and  the 
dynamic  free  surface  boundary  condition  contained  vortex  strength  in  the 
denominator.  A  check  was  therefore  performed  for  F  -  '  c-5  on  a  domain  sized 
five  ambient  wavelengths  in  both  depth  and  width.  The  vortices  were  located 
near  the  center  of  the  domain,  and  moved  little  during  the  checkout  time  scale. 

The  numerical  computations  performed  in  the  simulation  can  be  divided 
into  two  phases.  The  first  phase  concentrated  on  the  calculation  of  the  Cauchy 
integral  around  the  boundary.  In  order  to  test  the  accuracy  of  this  portion,  a 
deep  water  wave  velocity  potential  (<j>)  was  distributed  over  the  boundary 
elements  and  the  program  was  used  to  find  the  corresponding  stream  function 
(y).  This  test  simulated  the  first  step  in  the  final  program  format,  and  contained 
no  time  stepping  influence.  An  error  could  then  be  measured  by  comparing  the 
calculated  stream  function  and  the  exact  analytic  solution.  'P  z  maximum  error 
in  this  calculation  occurred  at  the  peaks  and  troughs  of  the  waves,  as  suspected, 
where  the  magnitude  of  the  second  derivative  of  the  potential  is  greatest.  A 
logarithmic  plot  of  this  error  versus  the  boundary  element  length  is  given  in  Fig 
6.  The  slope  of  this  plot  was  measured  as  2.16.  A  slope  of  2  a  esponds  to 
the  theoretical  error  limit  for  use  of  the  trapezoidal  rule  in  evaluation  of  the 
Cauchy  boundary  integral.  This  shows  that  trapezoidal  integration  error  is  the 
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dominant  numerical  approximation  in  the  calculations  and  could  be  reduced  by 
implementing  a  higher  order  integration  metho The  trapezoidal  method  of 
boundary  integration  was  retained  for  this  study  due  to  its  simplicity  however, 
and  a  sufficiency  fine  boundary  grid  size  was  chosen  to  evaluate  the  dominant 
character  of  the  vortex  interaction. 

The  second  phase  of  the  calculations  looked  at  the  time  stepping  procedure. 
The  velocity  potential  and  free  surface  elevation  were  integrated  forward  in  time 
using  the  free  surface  boundary  conditions  in  a  predictor/corrector  method.  In 
addition,  numerical  derivatives  of  the  stream  function,  velocity  potential  and 
surface  elevation  must  be  calculated  along  the  free  surface  for  inclusion  in  the 
boundary  condition  equations.  For  stability  reasons,  a  second  order  central 
difference  scheme  was  used  for  derivative  calculation  except  at  die  end  points, 
where  a  second  order  one  sided  difference  scheme  was  used.  Numerical  errors 
resulting  from  these  calculations  were  tested  by  combining  the  time  stepping  with 
the  Cauchy  boundary  integration  in  order  to  simulate  the  full  form  of  the 
program.  Again,  a  known  deep  water  wave  was  input  and  allowed  to  propagate 
one  wavelength  before  calculations  were  terminated.  The  calculated  free  surface 
profile  after  moving  one  wavelength  was  then  compared  to  the  exact  solution. 

The  errors  found  in  the  comparison  were  composed  of  both  dissipation  and 
dispersion  components.  These  can  be  seen  in  Fig  7  which  shows  the  surface 
elevation  solution  of  a  plane  travelling  wave  using  16  grid  points  per  wavelength. 
Over  the  time  scale  of  one  period  of  the  wave,  this  case  shows  a  4.8%  decrease 
in  the  mean  wave  amplitude,  and  a  1.7%  oscillation  in  amplitude  at  a  frequencv 
roughly  twice  the  wave  frequency.  This  error  is  the  manifestation  of  the  Cauchy 
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boundary  integration  errors  described  above,  the  numerical  differentiation  errors  in 
evaluation  of  the  potential  and  elevation  gradients  on  the  free  surface,  and  the 
time  integration  of  the  potential  and  elevation  using  the  boundary  condition 
equations.  Figure  8  shows  the  influence  of  varying  grid  density  on  the 
magnitude  of  these  errors.  A  resolution  of  64  grid  points  per  wavelength  was 
chosen  for  all  of  the  simulations  shown  in  this  report  except  the  ship  wake 
model.  The  domain  size  was  maximized  at  this  resolution  to  be  14  ambient 
wavelengths  wide.  This  domain  size  was  limited  by  the  available  core  memory 
of  the  DTRC  Cray  XMP  used  for  the  simulations.  Although  larger  memory  sizes 
would  allow  greater  size  and  accuracy  in  modelling  the  ship  wake,  the  current 
size  proved  adequate  for  the  purposes  of  this  study.  The  decay  over  one  wave 
period  at  this  density  was  0.25%,  and  the  oscillation  amplitude  was  approximately 
0.4%. 

3.2.3  Comparison  With  Nonlinear  Method 

Before  delving  into  the  complex  simulations  involving  ambient  waves, 
simulations  were  run  for  the  vortex  pair  without  ambient  waves.  The  purpose  of 
this  test  was  to  compare  simulation  output  with  independent  results  obtained  by 
Telste20  using  a  fully  nonlinear  free  surface  method.  From  these  comparisons, 
some  insight  into  the  limitations  of  the  linearized  method  could  be  determined, 
and  considered  in  interpretation  of  the  more  complicated  problem  involving 
ambient  wave  modification. 

Two  Froude  number  cases  were  used  for  comparison.  These  corresponded 
to  the  ’’low"  and  "moderate"  vortex  strength  cases  presented  by  Telste.  The 


"high"  Froude  number  case  was  not  simulated  because  nonlinear  vortex 
trajectories  in  this  case  were  significantly  influenced  by  the  free  surface  and  did 
not  closely  follow  the  path  expressly  used  in  the  linearized  method.  This 
limitation  in  the  linear  theory  may  be  partially  overcome  at  moderate  Froude 
number  values  by  utilizing  Cauchy’s  integral  theorem  to  compute  the  induced 
velocity  at  the  current  vortex  position  from  potentials  on  all  boundaries.  This 
additional  refinement  was  deemed  unnecessary  at  the  low  Froude  numbers  of 
concern  to  the  ship  wake  problem  of  interest  and  was  therefore  not  implemented. 
Insufficient  nonlinear  data  is  currently  available  to  accurately  evaluate  the  critical 
Froude  number  below  which  F  =  0  vortex  trajectories  are  sufficiently  accurate. 
The  results  show  however  that  Froude  numbers  below  Telste’s  "moderate"  value 
could  be  simulated  with  rigid  wall  vortex  paths.  In  the  cases  simulated,  the 
initial  separation  and  depth  of  the  pair  are  1.0  and  5.0,  respectively. 

The  "low"  Froude  number  condition  results  are  plotted  in  Fig  9.  Data  in 
this  figure  correspond  to  F  =  r/(gS3)w=  0.5  and  a  nondimensional  time  sufficient 
for  the  vortices  to  come  to  a  steady  translating  horizontal  speed  below  the 
surface.  The  x  axis  in  Fig  9  is  referenced  to  the  vortex  position  and 
nondimensionalized  by  the  initial  pair  separation  distance.  It  can  be  seen  that  the 
ALM  gives  very  good  agreement  with  the  nonlinear  method.  The  slight  under 
prediction  of  the  depth  of  the  surface  depression  is  most  likely  due  to  the 
increase  in  local  dynamic  pressure  in  this  region  caused  by  the  satisfaction  of  the 
surface  boundary  condition  on  y  =  0  in  the  linearized  case.  The  closer  proximity 
of  the  free  surface  to  the  vortex  in  the  nonlinear  method  results  in  a  larger 
induced  velocity  at  the  surface,  causing  a  lower  dynamic  pressure  and  further 
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surface  lowering  until  an  equilibrium  is  achieved.  This  figure  also  shows  the 
influence  of  neglecting  the  additional  terms  (basically  the  0.50 term)  in  the 
solution  obtained  using  the  CLM.  The  absence  of  induced  current  terms  prevents 
depression  of  the  surface  in  the  region  just  above  the  vortex  giving  an  upward 
hump  rather  than  a  dip  in  the  surface  above  the  vortex. 

A  comparison  of  surface  displacements  for  the  "moderate"  Froude  number 
is  shown  in  Fig  10.  The  data  correspond  to  F  =  2.236  and  t  =  32.  In  this  case, 
the  ALM  gives  qualitative  agreement  in  the  surface  profile  but  a  vertical  offset  in 
predicted  elevation.  This  offset  is  presumably  due  to  the  difference  in  position 
of  the  boundary  on  which  the  free  surface  conditions  are  satisfied  as  described  in 
the  previous  paragraph.  The  offset  error  at  this  Froude  number  is  not  considered 
detrimental  as  long  as  the  Froude  number  of  the  simulations  remain  below  this 
point.  This  condition  is  therefore  adopted  as  the  bound  in  simulation  conditions 
in  the  data  presented  in  this  study. 

Overall,  the  comparison  is  surprisingly  good  at  low  to  moderate  Froude 
numbers.  In  the  ambient  wave  simulations,  the  moderate  Froude  number 
described  above  was  not  exceeded  due  to  the  unknown  quality  of  surface 
deformation  results  in  this  range. 

3.2.4  Effects  of  Vortex  Induced  Current  on  Ambient  Waves 

The  maximum  induced  surface  velocity  from  a  single  vortex  and  its 
oppositely  signed  image  above  the  surface  can  be  written  as  v,  =  T/tcL,  where  L 
is  the  depth  of  the  vortex  from  the  surface.  The  ambient  wave  group  velocity  is 
defined  by  v4  =  0.5(gK/2n)m  for  linearized  deep  water  waves  (Newman11)-  The 


ratio  of  these  velocities  can  be  represented  using  the  previously  defined 
nondimensional  variables  in  the  following  way  (where  L  is  approximated  by 
S/2X): 

v,/vf  =  4  F/(S/X)  where  F  =  I7(2rtgX3)ia 

S  represents  1/2  the  depth  of  the  vortex  divided  by  the  ambient  wavelength  and 
its  use  here  is  a  precursor  to  the  geometric  conditions  of  the  vortex  pair  to  be 
studied  in  detail  in  the  next  chapter.  The  case  study  of  induced  current  effects 
involves  comparison  of  identical  flow  simulation  with  and  without  induced 
current  terms  in  the  free  surface  boundary  conditions.  Specifically,  the  0.5<J>02  + 
<j>0<l>i  terms  are  removed  from  the  dynamic  boundary  condition.  The  case  chosen 
for  this  comparison  contains  parameter  values  of  F  =  0.2  and  S/A.  =  1  and  is 
described  further  in  the  Froude  Number  Variation  section  later.  Figures  1 1  and 
12  show  the  corresponding  results  in  a  format  deserving  some  description.  The 
surface  profile  is  plotted  in  each  figure  for  a  series  of  times  corresponding  to  the 
period  of  the  undisturbed  ambient  wave.  Each  profile  is  offset  the  indicated 
amount  to  prevent  overlapping  and  facilitate  comparison  at  different  times. 
Horizontal  lines  are  drawn  for  each  profile  at  the  y  =  0  value  and  their  indicated 
spacing  provides  the  vertical  scale  for  wave  amplitude  measurement.  These 
figures  show  the  dramatic  effect  the  induced  current  has  on  changing  the  length 
of  the  ambient  waves  in  the  region  above  the  vortices.  Wavelength  here  is 
loosely  defined  as  the  distance  between  peaks  in  the  surface  elevation.  Without 
the  current,  the  maximum  stretched  wavelength  is  approximately  130+5%  of  the 
original  wavelength  at  t  =  13T  (T  =  period  of  ambient  wave).  The  induced 
current  case,  over  the  same  time  scale  produces  wavelengths  of  210+5%  the 
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original  wavelength.  The  shortest  waves  at  t  =  13T  are  80+5%  and  50+5%  of 
the  original  wavelength  in  cases  without  and  with  current  terms,  respectively. 

The  next  step  is  to  identify  how  the  parameters  of  a  certain  problem  (in 
this  case  the  vortex  pair)  influence  this  induced  current,  and  what  form  of  wave 
modifications  are  expected. 
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Figure  5  -  BIT  Domain  Description 
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Figure  10  -  Moderate  Froude  Number:  Comparison  of  Theories 
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Figure  11  -  Ambient  Wave  Modification:  Classical  Linearization 
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Figure  12  -  Ambient  Wave  Modification:  Alternate  Linearization 


CHAPTER  4  -  PARAMETRIC  STUDY 


4.1  NONDIMENSIONAL  PARAMETER  IDENTIFICATION 

The  incorporation  of  ambient  waves  into  the  problem  results  in  the  addition 
of  a  two  new  length  scales,  namely  the  length  and  amplitude  of  the  wave.  The 
previously  derived  boundary  conditions  need  coefficient  modification  to 
incorporate  these  additional  parameters.  Because  the  wavelength  is  directly 
related  to  the  linearized  wave  group  velocity  by  the  following  equation:  V,  = 
0.5(gXy2n)w,  the  group  velocity  could  be  held  constant  by  nondimensionalizing 
variables  based  on  the  ambient  wavelength  scale.  This  resulted  in  the  following 
nondimensionalization  scheme: 

x  =  x’A.  y  =  y'fk  t]  =  x\'fk  <J>  =  <j>7T  t  =  t’X2/T 
F=r/(gS3)in  P=A fk  Q=S  fk  R=S/D 

where  ’  indicates  dimensional  quantity,  and 
T  =  vortex  circulation 
S  =  inital  vortex  pair  separation 
D  =  initial  vortex  pair  depth 
X,  =  ambient  wavelength 
A  =  ambient  wave  amplitude 

Substituting  these  relations  into  the  previously  defined  ALM  boundary  conditions 
gives  the  following  results: 

kinematic  boundary  condition:  Tj,  =  <J>y  +  “n<J>0  -  <j)0Ti, 

”  yy  * 

dynamic  boundary  condition:  <|>t  =  -  0o0i  ■  fi/O^Q3) 
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The  corresponding  deepwater  wave  potential  using  the  above  nondimensional 
parameters  takes  the  form: 

=  P/(F  Q3fl  (2rc)lfl)  exp2**  sin  [2ttx  -  (27t)wt/(F  Q3/2)] 

The  vortex  potential  is  also  nondimensionalized  and  given  below  in  complex 
form: 

Po  =  <j>0  +  i\jr0  =  1/i  In  [(z-z^z-z^]  -  1/i  In  [(z-z*,)/(z-z*j)] 
where  the  character  represents  complex  conjugate.  This  collection  of 
equations  show  the  three  parameters,  F,  P  and  Q  as  intrinsic  coefficients.  F  is 
the  familiar  Froude  number  for  vortex  pairs  with  a  still  free  surface.  P,  or  AA, 
represents  the  steepness  ratio  of  the  ambient  waves.  Q,  or  S [k,  is  a  scaling 
parameter  between  the  vortex  spacing  and  the  ambient  wavelength.  R,  or  S/D, 
will  be  shown  later  to  represent  the  initial  position  of  the  vortices  on  their 
predefined  paths.  Since  S/D  and  SA  are  both  related  to  the  initial  geometric 
conditions  of  the  vortex  pair  their  influence  will  be  grouped  together.  Before 
proceeding  further,  these  geometric  parameters  and  the  associated  vortex  paths 
prescribed  will  be  described. 

As  a  consequence  of  the  linearization  as  discussed  in  the  previous  chapter, 
the  vortex  trajectories  follow  the  paths  of  F  =  0  solutions.  The  paths  prescribed 
for  this  rigid  free  surface  condition  were  first  described  by  Lamb7  as 

az(x2+y2)  =  4xY  where  "a"  is  an  arbitrary  constant.  (8) 

Paths  corresponding  to  this  equation  with  different  values  of  "a2"  are  plotted  in 
Fig  13.  It  should  be  noted  that  the  vortex  strength  is  excluded  from  this  path 
definition.  The  only  role  the  vortex  strength  plays  is  in  determining  the  speed  at 
which  the  vortices  move  along  thier  paths.  Since  the  vortices  must  follow  one  of 
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the  paths  shown,  their  positions  can  be  determined  by  knowing  their  start  point, 
their  strength,  and  the  elapsed  time.  This  figure  also  shows  that  the  final  depth 
of  the  vortices  is  equal  to  half  the  initial  separation  distance  if  the  initial  start 
point  was  rel?tively  deep  (small  S/D).  Cases  of  vortex  depth  less  than  half  of 
the  ambient  wavelength  were  avoided  due  to  the  inherent  errors  introduced  in 
vortex  paths  caused  by  wave  induced  velocities  at  these  shallow  depths. 

4.2  OUTPUT  FORMAT  DESCRIPTION 

The  output  of  the  transient  simulation  is  the  free  surface  profile  as  a 
function  of  the  nondimensional  time  variable.  The  quantity  of  output  generated 
is,  however,  somewhat  cumbersome  >  its  raw  form  due  to  the  large  number  of 
free  surface  grid  points  and  time  steps.  Three  dimensional  surfaces  composed  of 
the  2-dimensional  free  surface  profiles  stacked  in  rime  can  be  used  for  qualitative 
observation  of  the  data,  but  a  more  direct  2-dimensional  presentation  of  segments 
of  the  data  will  be  used  as  the  primary  presentation  form.  These  2-dimensional 
plots  will  show  the  surface  profile  at  intervals  of  time  corresponding  to  the 
penod  of  the  ambient  wave.  Each  surface  profile  is  offset  vertically  to 
differentiate  it  from  the  others  but  maintains  the  same  scale  for  comparison.  The 
sampling  of  the  output  at  periodic  intervals  results  in  very  little  qualitative  loss  in 
comprehension  of  the  flow  due  to  the  continuous  nature  of  the  surface  evolution 
as  verified  in  the  3-dimensional  surface  representations.  This  can  be  seen  in  Figs 
14  and  15  representing  a  typical  data  set.  Figure  14  shows  the  results  of  partial 
wave  reflection  on  the  right  side.  These  ambient  wave  rejections  form  a  pattern 
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similar  to  standing  waves.  Figure  14  also  shows  the  shortening  and  elongation 
process  as  a  bending  of  wave  crest  lines  in  time. 


43  FROUDE  NUMBER  VARIATION 

The  Froude  number  for  this  problem  can  be  defined  in  several  ways 
depending  on  one’s  choice  of  length  scale.  The  traditional  length  scale  chosen  in 
studies  without  ambient  waves  is  the  initial  vortex  separation  distance.  The 
resulting  form  of  the  Froude  number  is  F  =  I7(gSi)l/2.  The  examination  of 
Froude  number  influences  will  consist  of  flow  simulations  with  the  vortices 
started  at  a  relatively  deep  position  (S/D  =  1/3).  This  value  of  S/D  was  chosen 
to  minimize  surface  disturbances  associated  with  the  transients  of  the  ortices’ 
impulsive  start.  The  ambient  wave  modifications  would  therefore  h  associated 
predominately  with  the  velocity  field  generated  by  the  slowly  moving  vortices 
rather  than  the  initial  condition  wave  motion  generated  by  the  vortices.  S/D  = 

1/3  also  corresponds  well  with  the  experimental  studies  of  Willmarth  et.al.22 
described  earlier.  S/A.  was  set  equal  to  1  as  an  example  ratio.  This  value  of  S/A. 
allows  one  to  consider  the  Froude  number  to  be  based  on  the  length  scale  of  the 
ambient  waves  Fx  =  F  Q3'2  =  I7(gA,5)l/2  with  no  change  in  its  value.  The  resulting 
surface  profiles  for  F  ranging  between  0.1  and  0.6  are  shown  in  Figs  16  to  20. 
Since  the  time  scale  of  modification  is  much  slower  at  F  =  0.1,  the  profiles 
shown  in  Fig  16  coincide  with  2T  intervals  rather  than  the  IT  intervals  used  in 
Figs  17  to  20.  The  corresponding  vortex  paths  are  shown  for  comparison  in  Fig 
21.  Note  in  Fig  21  that  the  vortex  paths  of  each  condition  have  been  offset 
slightly  in  the  x  direction  to  lessen  marker  overlap.  In  actuality,  all  of  the  paths 
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are  the  same,  only  the  sample  points  along  the  path  differ.  The  surface 
modification  in  these  profiles  involves  a  stretching  of  waves  in  the  center  region 
above  the  vortices  at  later  times.  The  ambient  waves  also  appear  to  be  shortened 
somewhat  on  the  left  or  "windward"  side  due  to  the  adverse  directions  of  v,  and 
vr  The  wave  shortening  on  the  right  or  "leeward"  side  is  much  less  pronounced 
since  v,  and  vt  are  in  the  same  direction. 

The  most  obvious  effect  of  varying  the  F  in  these  cases  is  the  change  in 
time  scale  for  the  surface  disturbance  to  develop.  The  occurrence  of  ambient 
wave  modification  at  earlier  times  at  higher  F  is  due  primarily  to  the  increase  in 
relative  speed  of  the  vortex  pair  moving  towards  the  surface.  To  minimize  this 
attribute,  surface  wave  profiles  can  be  compared  between  different  F  cases  when 
the  vortex  pair  are  at  common  posldons.  Two  common  positions  are  labeled  "A” 
and  "B"  in  Fig  21.  The  corresponding  profiles  from  each  F  arc  collected  and 
shown  in  expanded  vertical  scale  in  Figs  22  and  23.  Case  "A"  corresponds  to  a 
moderately  deep  vortex  pair.  Only  the  higher  F  cases  corresponding  to  relatively 
stronger  vortices  have  an  influence  at  this  depth,  verifying  the  notion  that 
stronger  vortices  will  influence  the  surface  from  greater  depths. 

A  more  subtle  and  interesting  difference  is  observed  as  die  vortex  pair 
draws  near  the  surface  at  position  "B”.  Here  one  notices  the  lengthening  of 
residual  ambient  waves  in  the  central  region  above  the  vortices.  This  increase  in 
residual  wave  length  is  presumably  caused  by  the  higher  tangential  velocity 
gradients  at  the  surface  in  cases  where  F  is  increased.  The  width  of  the 
disturbed  region  above  the  vortex  pair  does  not  depend  heavily  on  the  value  of 
F,.  Figure  24  shows  the  vortex  induced  current  for  condition  B  without  waves. 


47 


For  threshold  values  of  vt/vt  above  0.25,  the  trend  in  width  of  effected  region 
increases  at  least  linearly  with  F,.  The  left  side  position  of  wave  compression  in 
Fig  23  suggests  that  the  real  growth  of  the  effected  region  is  less  than  linear. 
This  discrepancy  is  most  likely  due  to  the  time  lag  in  free  surface  response 
which  increases  as  the  vortex  speed  is  increased. 

An  additional  note  on  the  variation  of  F  is  the  transition  occurring  around 
values  of  0.5  of  the  surface  disturbance  directly  above  the  vortex  pair.  At  values 
above  0.5,  the  creation  of  a  center  "hump"  appears  similar  to  center  humps 
calculated  by  Telste20  at  higher  Froude  numbers  without  ambient  waves.  It 
should  be  noted  that  in  this  case  the  vortices  have  not  deviated  from  their  rigid 
wall  paths.  This  indicates  a  transition  Froude  number  value  of  0.5  above  which 
the  vortex  pair  begins  dominating  the  surface  profile. 

4.4  INITIAL  GEOMETRIC  CONDITIONS 

The  initial  position  of  the  vortex  pair  is  equally  important  to  its  strength  in 
defining  the  form  of  ambient  wave  modification.  As  noted  previously,  the  ratio 
of  induced  surface  velocity  to  ambient  wave  group  velocity  is  highly  dependant 
on  the  depth  of  the  vortex.  Time  scales  for  the  vortex  pair  to  approach  the 
surface  from  large  depths  is  also  obviously  dependant  on  the  initial  depth  and 
separation  distance.  The  paths  of  the  vortices  in  all  of  the  cases  studied  were 
defined  by  the  rigid  surface  condition  as  explained  in  Chapter  I.  This  family  of 
paths  described  by  equation  8  and  shown  in  Fig  13  show  that  the  position  of 
vortex  pair  can  be  identified  by  two  parameters.  The  parameter  S/D  specifies  the 
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relative  position  of  the  vortex  on  any  of  the  family  of  paths.  The  particular  path 
desired  may  then  be  specified  by  designating  a  value  for  either  S  or  D. 

4.4.1  S/D  Variation 

The  influence  of  starting  at  various  points  on  a  given  path  will  be 
examined  first.  The  controlling  parameter  in  this  case  is  S/D,  which  determines 
the  initial  point  on  the  path.  Three  values  of  S/D  are  examined  corresponding  to 
deep  (S/D  =  0.5),  turning  (S/D  =  2),  and  shallow  (S/D  =  6).  As  shown  earlier, 
the  direction  of  motion  of  the  vortex  pair  is  quite  different  in  each  of  these 
cases.  The  value  of  F*  used  in  these  cases  was  set  to  0.125,  and  the  common 
path  chosen  corresponded  to  a  final  submergence  depth  of  the  vortices  of  Dw  = 
X.  Figures  25  to  27  show  the  resulting  profiles  and  Fig  28  shows  the 
corresponding  trajectories  of  the  right  side  vortex. 

4.4.2  S fk  Variation 

As  mentioned  in  the  S/D  section  above,  initial  vortex  pair  location  is 
important  in  examining  the  influences  on  ambient  surface  waves.  In  addition  to 
specifying  the  portion  of  the  path  that  the  pair  will  be  started  on,  a  particular 
path  must  be  chosen.  The  parameter  S fk  is  used  as  a  path  selection  parameter 
for  its  role  in  scaling  the  vortex  path  geometry  to  the  ambient  wavelength. 

Three  values  of  S (k  were  simulated  with  a  constant  value  of  S/D  =  6.  S/X 
values  were  set  to  3.0,  4.5,  and  6.0.  The  upper  limit  on  S/X  is  caused  by  the 
size  limitation  of  the  computational  domain,  and  the  influence  of  partial  wave 
reflections  on  the  right  hand  boundary.  Fx  varied  between  0.06  and  0.125  to 
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preserve  a  \  J\t  at  the  surface  of  0.5.  Surface  profiles  for  these  simulations  are 
shown  in  Figs  29  to  31.  The  corresponding  vortex  trajectories  are  given  in  Fig 
32. 


4.4.3  Comments 

The  profiles  shown  in  this  study  of  geometric  conditions  on  the  vortex  pair 
indicate  the  great  complexity  of  the  flow  described  by  the  equations  outlined 
previously.  The  detailed  description  of  the  flow  field  and  its  dependence  on 
these  parameters  is  therefore  unpractical  at  this  stage  of  development.  Several 
very  important  general  characteristics  of  the  ambient  wave  modification  can  be 
described,  however,  providing  greater  understanding  of  the  important  physical 
mechanisms  governing  the  interaction. 

The  first  general  characteristic  found  was  the  dominance  of  the 
instantaneous  depth  of  the  vortex  pair.  This  is  seen  primarily  in  the  difference  in 
surface  disturbance  caused  by  the  vortex  pair  initialized  at  different  points  on  the 
same  path  (varying  S/D).  In  the  case  of  S/D  =  0.5,  the  surface  is  modified  only 
by  slight  shifting  in  the  waves  caused  by  the  small  vortex  induced  surface 
current.  This  wave  alteration  is  also  noticed  to  be  almost  symmetric  about  the 
line  of  symmetry  for  the  vortex  paths.  This  indicates  little  net  exchange  of 
energy  between  the  vortices  and  the  waves  as  they  pass  through.  As  the  depth  is 
decreased,  the  vortex  induced  surface  current  grows  to  the  value  of  the  group 
velocity  of  the  ambient  waves,  and  essentially  stops  the  transfer  of  ambient  wave 
energy  from  left  to  right.  Two  important  processes  occur  which  deserve  further 
attention.  First,  the  ambient  waves  are  shortened  as  they  enter  the  region  of 


50 


influence  of  the  vortex  pair  by  the  apposing  induced  current.  The  resultant 
decrease  in  ambient  wavelength  is  accompanied  by  a  decrease  in  the  absolute 
group  velocity  of  the  waves.  The  result  of  this  process,  is  that  ambient  waves  of 
higher  speed  may  be  halted  by  vortices  of  smaller  magnitude  than  that  predicted 
from  vi/\t  =  1.  The  determination  of  either  a  critical  vortex  strength  or  ambient 
wave  length  for  wave  stoppage  is  frustrated  by  the  influence  of  time  scales 
associated  with  the  motion  of  the  vortices  and  their  depth.  In  addition,  the 
absolute  identification  of  the  wave  blockage  condition  threshold  is  not  currently 
feasible. 

A  second  important  process  discovered  in  the  simulations  is  the  dynamics 
of  waves  caught  in  the  central  influence  region.  These  surviving  center  waves 
may  be  left  over  from  the  initial  ambient  wave  profile  in  this  region  in  the 
shorter  time  scales.  These  left  over  waves  are  found  in  greater  numbers  as  Sfk 
is  increased,  as  shown  in  Figs  29  to  31.  The  cases  shown  in  these  figures  show 
a  quick  generation  of  wave  blocking  action  due  to  the  relatively  large  value  of  vt 
/vr  Wave  energy  may  also  propagate  into  the  center  region  when  the  vortex 
pairs  are  initialized  at  low  values  of  S/D.  This  energy  propagation  is  facilitated 
by  the  cancelling  influence  of  the  induced  tangential  velocity  from  each  of  the 
counter  rotating  submerged  vortices  at  low  values  of  S/D.  The  evolution  of  these 
center  waves  is  particularly  complex  in  cases  of  moderate  S/D.  Initial  waves  in 
this  region  are  typically  stretched  by  the  surface  current  gradient  The  resulting 
increase  in  wavelength  gives  rise  to  an  increase  in  phase  and  group  velocity. 
Because  the  gravity  wave  dispersion  relation  specifies  that  longer  waves  move  at 
higher  velocity,  the  central  region  may  be  evacuated  over  a  shorter  time  than  that 


51 


predicted  by  ambient  wave  group  velocity.  Simulation  results  were  not  carried 
far  enough  in  time  to  show  complete  smoothing  of  the  central  region  due  to  the 
influence  of  partial  wave  reflection  from  the  right  side  boundary.  Larger  surface 
grids  and  greater  computer  time  will  therefore  be  necessary  to  validate  this 
hypothesis. 

4.5  AA  VARIATION 

The  parameter  AA  found  in  the  derivation  of  the  ALM  represents  the 
influence  of  ambient  wave  steepness  on  the  interaction  process  with  submerged 
vortices.  To  examine  its  influence  a  case  of  AA  =  1/10  was  simulated  and 
compared  with  the  same  case  at  the  baseline  AA  =  1/20.  The  resulting  profiles 
from  these  two  cases  are  shown  in  Fig  33.  The  conditions  for  the  comparison 
were  Fx  =  0.125,  S/D  =  2,  and  SA  =  2.  No  apparent  changes  in  the  results 
other  than  a  constant  scale  factor  are  present.  Additional  analysis  (Fig  34)  shows 
the  result  of  multiplication  of  the  AA  =  1/20  elevation  results  by  2  and 
subtracting  the  AA  =  1/10  results.  From  this  plot  one  notices  the  weak  coupling 
of  this  parameter  with  the  long  wavelength  character  of  the  waves  generated  by 
the  initial  conditions  of  the  vortex.  This  is  verified  by  almost  identical  results 
between  Fig  34  and  simulations  with  AA  =  0.  It  was  concluded  from  this  that 
the  wave  steepness  parameter  plays  a  minor  role  in  shaping  the  wave 
modification  region,  and  provides  a  roughly  linear  vertical  amplification  of 
simulation  profiles. 
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4.6  PHASE  VARIATION 


Before  pursuing  spectral  analysis  of  the  simulation  results,  the  influence  of 
ambient  wave  phase  must  be  scrutinized.  A  sample  case  is  simulated  with 
multiple  phase  ambient  waves  and  the  resultant  surface  profile  shown  in  Fig  35. 
This  figure  shows  no  significant  variation  in  wave  modification  as  the  phase  is 
changed.  In  addition,  the  spectral  analysis  described  in  the  next  section  was 
applied  to  each  of  these  cases  with  no  significant  variation  associated  with  the 
ambient  wave  phase.  The  influence  of  ambient  wave  phase  will  therefore  be 
considered  secondary  to  the  other  parameters  of  the  problem  and  will  be 
neglected  in  the  spectral  analysis  to  follow. 


\ 
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Figure  14  -  3-D  Free  Surface  Representation 
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Figure  15  -  Discrete  Profile  Surface  Representation 
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Figure  17  -  F  Variation:  F=0.2,  S/D=0.33,  S/X=l 
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Figure  18  -  F  Variation:  F=0.4,  S/D=0.33,  S/X=l 
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Figure  20  -  F  Variation 


:  F=0.6,  S/D=0-33,  Sfk-) 
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Figure  22  -  F  Variation:  Position  "A"  Surface  Profile  Comparison 
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Figure  23  -  F  Variation:  Position  "B"  Surface  Profile  Comparison 


62 


X/x 


Figure  24  -  F  Variation:  Induced  Surface  Current  at  Position  "B" 
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Variation'.  S/D=2.0,  F=0.125,  S/X=2.5 
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Figure  29  -  Sfk  Variation:  S/X=3.0,  F=0.06,  S/D=6.0 
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Figure  31  <  S/A.  Variation:  SA=6.0,  F=0.125,  S/D=6.0 
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Figure  33  Ambient  Wave  Amplitude  Variation 
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CHAPTER  5  -  SPECTRAL  ANALYSIS 


In  order  to  make  the  previous  findings  useful  in  the  study  of  real  ocean 
flows,  a  method  of  analysis  must  be  developed  which  illuminates  the  impact  of 
wave  modification  on  radar  images.  In  Chapter  1,  it  was  shown  that  the 
reflected  radar  intensity  is  proportional  to  the  amount  of  Bragg  waves  present  in 
the  region  being  illuminated  by  the  radar.  This  intensity  is  usually  expressed  as 
the  backscattering  cross  section  per  unit  surface  area,  o0.  This  cross  section  is 
derived  by  Valenzuela  (1978)  in  the  well  known  Bragg  scattering  expression: 
o0  =  4nkI4cos40  F(0)  M(2kTsin0) 

where:  k,  =  radar  wavenumber  and  0  =  vertical  incidence  angle  of  radar.  Here 
F(0)  is  a  scattering  coefficient  dependant  on  the  radar  polarization  and  incidence 
angle,  and  M()  is  the  component  of  the  wave  spectra  corresponding  to  the  Bragg 
condition.  Remember,  from  the  introduction,  that  the  Bragg  condition  is  satisfied 
when  k*  (surface  wavenumber  component  in  the  look  direction  of  the  radar)  is 
equal  to  the  argument  of  M()  in  the  equation  above.  If  0  and  k,  are  known 
constants,  then  o0  is  proportional  to  M().  This  property  indicates  that  a  relative 
variation  in  radar  backscatter  over  different  sections  of  the  water  surface  can  be 
determined  by  the  variations  in  the  Bragg  wave  spectral  component  among  the 
respective  surface  sections.  Estimates  of  the  relative  radar  intensity  image  can 
therefore  be  compiled  by  performing  Finite  Fourier  Transforms  of  the  surface 
elevation  over  regions  corresponding  in  size  to  the  resolution  cell  size  of  the 
radar.  This  will  be  performed  here  by  sweeping  a  window  across  the  simulated 
surface  data  and  calculating  the  Finite  Fourier  Transform  at  various  window 
positions.  The  Fourier  Transform  component  corresponding  to  the  Bragg  waves 


in  each  window  position  is  then  assigned  to  the  midpoint  location  of  the  window. 
The  resulting  distribution  of  Bragg  wave  components  over  the  surface  is  then 
contour  plotted  for  comparison  with  SAR  imagery. 

5.1  METHOD 

The  Finite  Fourier  Transform  used  to  determine  the  Bragg  wave  spectral 
content  in  a  data  window  is  given  by: 

F(k,L)  =  I  f(x)  el2*kxdx 

-L/2-' 

where:  L  is  the  window  width.  The  application  of  this  equation  to  discrete  data 
can  be  performed  using  an  algorithm  known  as  the  Fast  Fourier  Transform,  or 
FFT.  This  common  method  is  described  in  detail  by  Bendat  and  Piersol.1  The 
FFT  output  coefficients  are  complex  in  form  providing  phase  information  for 
each  wavenumber  component.  Since  this  phase  was  shown  to  be  unimportant  in 
the  previous  chapter,  the  absolute  magnitude  of  each  wavenumber  component  is 
formed  from  the  modulus  of  the  real  and  imaginary  parts: 

M(k)=[FFT^(k)  +  FFTaM(J(k)]1'2 
where:  FFT(k)  is  an  output  coefficient. 

Before  applying  the  scanning  FFT  window  to  simulation  results,  several 
comments  must  be  made  here  to  avoid  misinterpretation  of  the  method.  The 
regions  used  in  computation  of  the  spectra  will  be  overlapped  here  because  of  the 
finite  extent  of  the  surface  domain  and  the  desire  to  have  moderate  sized  surface 
patches.  If  the  patches  are  allowed  to  be  too  small,  the  resolution  of  low 
wavenumber  components  decreases,  resulting  in  undistinguishable  movements  of 
wave  energy  among  the  wavenumbers  of  concern.  This  overlapping  is  not 
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representative  of  radar  imaging,  and  therefore  the  results  must  not  be  considered 
direct  simulation  of  radar  operation.  The  general  trends  in  the  resultant 
component  levels  in  the  Bragg  wavenumber  range  are  believed  to  be  somewhat 
representative  of  results  obtained  with  radar  processing. 

5.2  APPLICATION  TO  SHIP  WAKE 

The  case  examined  using  this  technique  is  a  simulation  approximating  the 
wake  velocities  from  a  twin  screw  destroyer  calculated  by  Sween  (1987).  It 
should  be  noted  that  in  simulating  the  full  ship  case,  the  grid  resolution  was 
lowered  considerably  by  constraints  in  available  computer  memory.  The  grid 
density  is  32  points  per  ambient  wavelength.  The  ambient  wavelength  is  set  to 
25  cm  and  corresponds  to  the  Bragg  scattering  wavelength  for  an  "L  band"  SAR 
operating  at  an  incidence  angle  of  30  degrees.  The  initial  vortex  separation  and 
depth  are  5  m  and  2  m,  respectively.  Gravitational  acceleration  is  equal  to  9.8 
m/s2,  and  the  ambient  wave  amplitude  is  set  to  2.5  cm.  The  vortex  circulation 
strength,  T,  is  chosen  to  give  approximately  the  same  maximum  surface  current 
reported  by  Sween  (1987)  for  the  given  depth  of  submergence.  The  resulting  T 
value  was  1.0.  The  simulation  domain  extended  to  +  6.25  m  in  x  direction  with 
a  depth  of  12.5  m.  The  simulation  was  carried  out  for  8  periods  of  the  ambient 
wave  (3.2  s)  and  the  resultant  profiles  at  each  half  period  are  shown  in  Fig  36. 
The  corresponding  downstream  distance  spanned  by  the  simulation  is  25  m  for  a 
ship  speed  of  15  knots. 

The  FFT  window  size  was  chosen  to  be  2  m  wide,  which  is  typical  of 
SAR  resolution  cell  dimensions.  The  resultant  spectral  resolution  is  3.14/cm  in 
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wavenumber  giving  components  at  wavelengths  of  22  cm,  25  cm,  and  28  cm. 

The  consecutive  FFT  windows  were  overlapped  in  the  analysis  by  3  X.  This 
overlap  gives  denser  data  for  examining  the  variations  in  Bragg  component, 
though  only  non-overlapping  sections  would  be  directly  comparable  to  typical 
radar  imaging  techniques.  The  Bragg  component  of  each  spectra  is  assigned  to 
the  location  of  the  FFT  window  center.  The  resulting  contour  plot  of  the 
magnitude  of  the  Bragg  component  versus  space  and  time  is  given  in  Fig  37. 

This  plot  indicates  a  reduction  in  Bragg  wave  component  in  the  center  region  as 
time  progresses.  In  order  to  understand  the  energy  flow  path  in  this  region.  Fig 
38  is  constructed.  This  figure  shows  the  progression  of  the  spectra  as  the 
window  is  traversed  across  the  contour  plot  at  the  frame  position  200.  The  thick 
line  indicates  the  wavenumber  associated  with  the  Bragg  wavelength.  Figure  38 
illuminates  the  shift  in  wave  energy  to  lower  wavenumbers  in  the  center  region. 
Since  these  wavenumber  components  do  not  satisfy  the  Bragg  scattering  criterion, 
the  resultant  received  energy  at  the  radar  is  decreased.  This  phenomenon  of 
wavenumber  shifting  due  to  induced  surface  currents  appears  to  produce  a  weaker 
radar  return  signal  in  the  region  between  the  vortices.  In  the  ocean,  the  wave 
field  can  be  described  by  the  superposition  of  many  wavelength  waves.  Typical 
spectra  from  the  sea  (Newman11)  indicate  a  decreasing  level  of  energy  as 
wavenumber  is  increased  in  the  range  of  the  L  Band  Bragg  scattering  waves. 

The  implication  of  this  is  that  though  higher  wavenumber  (shorter)  waves  will  be 
stretched  into  the  Bragg  wavelength,  thier  original  energy  is  lower.  Thus  one 
would  expect  that  a  larger  amount  of  energy  was  shifted  out  of  the  Bragg 
sensitivity  window  than  were  being  shifted  in  The  inclusion  of  these  higher 
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wavenumbers  would  however,  decrease  the  magnitude  of  variation  in  M()  shown 
in  Fig  37.  The  importance  of  this  action  is  that  it  produces  the  proper  trend  in 
the  surface  response  in  the  wake  region  of  a  ship  and  may  be  a  significant 
contributor  to  the  overall  signature  production  process  for  the  dark  centerline 
wake. 
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Figure  37  -  Ambient  Wavenumber  Fourier  Component  Space/Time  Distribution 


81 


Time  (seconds) 


Content  in 


CHAPTER  6  -  SUMMARY 


The  modification  of  surface  waves  by  submerged  vortices  has  been  explored 
here  in  a  straight  forward  progression  from  analytical  examination  through 
numerical  simulation  of  a  simple  vortex  ambient  wave  model.  The  motivation 
for  study  of  this  problem  is  the  determination  of  ship  wake  impacts  on  synthetic 
aperture  radar  images  showing  dark  centerline  regions  behind  ships  moving  at 
sea.  This  process  has  lead  to  several  discoveries  in  the  nature  of  the  interaction 
process  between  vortices  and  ambient  waves. 

In  the  initial  analytical  study,  a  new  form  of  linearized  equations  was 
derived  for  the  kinematic  and  dynamic  free  surface  boundary  conditions  involving 
submerged  vortices.  The  additional  terms  resulting  from  this  analysis  were 
shown  to  be  of  fundamental  importance  in  the  study  of  surface  flows  involving 
low  Froude  number  vortical  disturbances.  Comparisons  of  the  Alternate 
Linearization  Method  (ALM)  with  fully  nonlinear  solutions  show  very  good 
agreement.  Simulation  results  using  the  ALM  indicate  that  the  physical 
mechanism  responsible  for  significant  surface  disturbances  due  to  vortices  is  the 
vortex  induced  surface  current  Example  simulations  from  this  program  show 
that  neglect  of  boundary  condition  terms  representing  the  induced  current  give 
much  smaller  modification  of  surface  waves  than  equivalent  conditions  with  the 
induced  surface  current  included. 

The  simulation  program  has  also  been  used  to  examine  the  specific  flow 
configuration  of  a  vortex  pair  impulsively  started  below  a  free  surface  containing 
a  sinusoidal  ambient  wave  train.  This  examination  determines  which  parameters 
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are  significant  in  determining  the  form  of  the  ambient  wave  modification.  It  is 
found  that  the  Froude  number  of  the  vortices  and  their  position  relative  both  to 
each  other  and  the  free  surface  play  the  dominant  toles  in  modifying  ambient 
waves.  The  ambient  wave  steepness  and  initial  relative  phase  with  respect  to  the 
vortices  have  only  secondary  influence  in  the  surface  modification. 

The  simulation  is  then  extended  to  represent  the  simple  twin  vortex  wake 
model  of  a  generic  ship  wake.  The  results  of  this  simulation  are  analyzed  using 
a  scanning  Fast  Fourier  Transform  window  on  the  surface  to  obtain  spacial 
distributions  of  the  spectral  content  of  the  surface.  These  spacial  distributions  are 
filtered  to  extract  the  Bragg  scattering  wave  component  and  contour  plotted.  The 
resulting  image  shows  that  the  shift  in  Bragg  wave  energy  to  longer  wavelengths 
provides  a  mechanism  for  the  observance  of  darker  regions  in  centerline  wake 
regions  of  ships  moving  in  radar  images. 


REFERENCES 


1.  Bendat,  J.S.,  and  A.G.  Piersol,  Random  Data:  Analysis  and  Measurement 
Proceedures,  Wiley  Interscience,  New  York,  pp.  252-253  and  pp.  299-309, 
(1971). 

2.  Cooper,  AC.  "Interactions  Between  Ocean  Surface  Waves  and  Currents,”  NRL 
Mem.  Rep.  5755,  (1986). 

3.  Fish,  S.,  "Ambient  Free  Surface  Wave  Modification  by  a  Submerged  Vortex 
Pair,"  Ph.D.  Thesis,  Univ.  of  MD,  College  Park,  (1989). 

4.  Fish,  S.  and  J.N.  Blanton,  "Analysis  of  Propeller  Wake  Flow  Visualization 
Near  a  Free  Surface,"  DTRC/SHD  report  1268-03,  (1988). 

5.  Fry,  D.  and  T.T.  Huang,  "Turbulent  Ship  Wake  and  Ambient  Wave 
Interaction,"  IRIIED  Annual  Report,  DTRC  Jan.  1989,  pp  27-32. 

6.  Griffin,  O.M.  and  G.A.  Keramidas,  T.F.  Sween,  and  H.T.  Wang,  "Ocean  and 
Ship  Wave  Modification  by  a  Surface  Wake  Flow  Pattern,”  NRL  Mem.  Rep. 

6094,  (1988). 

7.  Lamb,  H.  Hydrodynamics,  6th  Ed.,  Cambridge  University  Press,  N.Y.  pp. 
223-224,  (1945) 

8.  Lindenmuth,  B.,  and  D.  Fry,  "Viscous  Macro-Wake  Behind  a  Twin  Screw 
High  Speed  Surface  Ship,"  DTRC  SHD  report  127301,  (1990). 

9.  Lyden,  J.,  D.  Lyzenga,  R.  Shuchman,  and  E.  Kasischke,  "Analysis  of  Narrow 
Ship  Wakes  in  Georgia  Strait  SAR  Data,"  ERIM  report  155900-20-T,  Ann  Arbor, 
MI.,  (1985). 

10.  Marcus,  D.L,  "The  Interaction  Between  a  Pair  of  Counter-Rotating  Vortices 
and  a  Free  Surface,"  Ph.D.  Thesis,  Dept,  of  Mech.  Eng.,  Univ.  of  Calif,  at 
Berkeley,  (1988). 

11.  Newman,  J.N.,  Marine  Hydrodynamics,  MIT  Press,  Cambridge,  MA.  pp. 
240-259,  (1977). 

12.  Ohring,  S.  and  H.J.  Lugt,  "Two  Counter-Rotating  Vortices  Approaching  a 
Free  Surface  in  a  Viscous  Fluid,"  DTRC  report  89/013,  (1989). 

13.  Peltzer,  R.,  and  W.  Garrett,  and  P.  Smith,  "A  Remote  Sensing  Study  of  a 
Surface  Ship  Wake,"  Int.  J.  Remote  Sensing  Vol  8,  no.  5,  pp.  689-704,  (1987). 


85 


14.  Sarpkaya,  T.  and  D.O.  Henderson  Jr.  "Free  Surface  Scars  and  Striations  Due 
to  Trailing  Vortices  Generated  by  a  Submerged  Lifting  Surface,"  paper  85-0445 
at  AIAA  23rd  Aero.  Sci.  Meeting,  Jan.  14-17,  Reno,  Nev.  (1985). 

15.  Sarpkaya,  T.  "Trailing- Vortex  Wakes  on  the  Free  Surface,"  16th  Symposium 
on  Naval  Hydrodynamics,  National  Academy  Press,  pp.  38-50,  (1986). 

16.  Sarpkaya,  T.,  Elinsky  II,  J.,  and  R.  E.  Leeker  Jr.  "Wake  of  a  Vortex  Pair  on 
the  Free  Surface,"  17th  Symposium  on  Naval  Hydrodynamics,  Aug.  29-Sep.  2,  The 
Hague,  The  Netherlands.,  (1988). 

17.  Simmen,  J.A.  and  P.G.  Saffman  "Steady  Deep-Water  Waves  on  a  Linear 
Shear  Current,"  Stud.  Appl.  Maths.  Vol  73,  pp.  35-57,  (1985). 

18.  Teles  Da  Silva,  A.F.  and  D.H.  Peregrine,  "Steep  Steady  Surface  Waves  on 
Water  of  Finite  Depth  and  Constant  Vorticity,"  J.  Fluid  Mech.  Vol  195,  pp.28 1- 
302.  (1988). 

19.  Tryggvason,  G.,  "Deformation  of  a  Free  Surface  as  a  Result  of  Vortical 
Hows,"  Phys.  Fluids  31  (5),  pp.  955,  (1988). 

20.  Telste,  J.G.,  "Potential  How  About  Two  Counter-Rotating  Vortices 
Approaching  a  Free  Surface,"  J.  Fluid  Mech.  Vol  201,  (1989). 

21.  Wehausen,  J.V.,  and  E.V.  Laitone,  "Surface  Waves,"  Handbuch  drr  Physik 
IX,  Springer- Verlag,  (1960). 

22.  Willmarth,  W.W.  and  G.  Tryggvason  and  A.  Hirsa  and  D.  Yu,  "Vortex  Pair 
Generation  and  Interaction  with  a  Free  Surface,"  Phys.  Fluids  A,  Vol  1,  Feb. 

1989,  pp.  170-172. 

23.  Wright,  J.W.  "Backscattering  from  Capillary  Waves  with  Application  to  Sea 
Clutter,"  IEEE  Trans.  AP-14,  pp.  749-754,  (1966). 


86 


INITIAL  DISTRIBUTION 


Copies 


Copies 


1  CHONR/Code  432  Rood 

1  NAVPGSCOL  Sarpkaya 

2  NRL 

1  Lib 
1  Griffin 

2  NAVSEA 

1  SEA05R 
1  SEA55W3 

1  NAVFACENGCOM 

1  NAVOCEANO/lib 

1  NADC 

1  NWC 

1  NOSC 

1  CElVCode  L31 

1  NSWC/White  Oak/Lib 

1  NSW  C/Dahl  gren/Lib 

1  NUSC  NPT 

12  DTIC 

2  AFFDLVFYS 

1  Dale  Cooley 
1  S.J.  Pollock 

1  COGARD 

1  COM  (E),  STA  5-2 

1  Long  Island  Univ 

Grad  Dept  of  Marine  Science 
David  Price 

1  /\FFDL/FDDS/J.  Olsen 


1  LC/Sci  &  Tech  Div 

1  MARAD/Adv  Ship  Prog  Office 

1  MMA/Tech  Lib 

1  NASA  AMES  RESEARCH  CTR 

R. T.  Medan,  Ms  221-2 

1  NASA  LANGLEY  RESEARCH 

D.  Bushnell 

1  NASA/Sci  &  Tech  Info  Facility 

1  NSF/Eng  Div 

1  Univ.  of  California,  Berkeley 

College  of  Eng,  NA  Dept 

S.  Berger 

1  CA  Inst  of  Tech 

A.J.  Acosta 

1  Cornell  University 

Grad  School  of  Aero  Eng 

1  Florida  Atlantic  Univ 

Ocean  Eng  Lib 

1  Univ  of  Hawaii/Dr.  Bretschneider 

1  Univ  of  Illinois/Coll  of  Eng 

J.M.  Robertson 
Theoretical  &  Applied  Mech 

1  State  Univ  of  Iowa 

Iowa  Inst  of  Hyd  Research 
L.  Landweber 

1  University  of  Maryland 

Eng  Lib 

1  VPI 

J.  Schetz 


87 


Copies 


Copies 


2 


4 


1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

2 


MIT/Dept  of  Ocean  Eng  1 

1  J.  Milgram 
1  J.N.  Newman 

1 

Univ  of  Michigan/Dept/NAME 
1  Lib 

1  R.  Beck  1 

1  Willmarth 
1  Trygvasson 

Penn  State  Univ  1 

Applied  Research  Lab 


SAl/Annapolis 
N.  Salvesen 


Southwest  Research  Inst 
Applied  Mech  Review 

Stanford  University 
Department  of  Div  Eng 
Dept  of  Aero  and  Astro 
M.  Van  Dyke 


Stanford  Research  Inst/Lib 

Webb  Institute 
Lib 


Copies 

1 

1 

1 

1 

1 

1 

1 

1 

5 

1 


Woods  Hole,  Ocean  Eng  Dept  1 

1 

SNAME  1 

Bolt,  Baranek  &  Newman,  MA  1 

1 

Boeing  Company/Aerospace  Group  jq 
R.R.  Barberr 


CALSPAN,  Inc.  Applied  Mech  Dept 
Bow  Research,  Inc. 


Eastern  Research  Group 

General  Dynamics  Corp 
1  Convair  Aerospace  Div 

A.M.  Cunningham,  Jr.  Ms  2851 
1  Electric  Boat  Div 
V.T.  Boatwright,  Jr. 


Gibbs  &  Cox,  Inc. 

Tech  Info  Control  Selection 

Tracor  Hydronautics,  Inc. 

A.  Goodman 

Lockheed  Aircraft  Corp 
Lockheed  Missiles  &  Space 
R.L.  Waid 

McDonnell-Douglas  Corp/Douglas 
Aircraft  Company  -  Lib 

CENTER  DISTRIBUTION 

Code  Name 
012 
12 

14  Ken  Nicolas 

15 
152 
154 

1542  Noblesse 

1543 

1543  Fish 

1544  Reed 
156 

1802  Lugt 
1843 

342.1  TIC  (C) 

342.2  TIC  (A) 

3432  Reports  Control 


88 


